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We discuss in detail the photon spectrum of radiative T(1S) decays taking into account a number 
of results that have recently appeared in the literature. In particular, we show how to consistently 
combine expressions which are valid in the upper end-point region, where NRQCD factorization 
breaks down, with those of the central region, where NRQCD factorization holds. An excellent 
description of data is achieved, but theoretical errors are large. 
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I. INTRODUCTION 

Semi-inclusive radiative decays of heavy quarkonium systems (see [l| for a review) to light hadrons have been a 

n 
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' subject of investigation since the early days of QCD 0,13. It was thought for some time that a reliable extraction 
of a s was possible from the photon spectrum normalized, for instance, to the decay into muon pairs. However, the 
upper end-point region of the spectrum (namely z —* 1, z being the fraction of the maximum energy the photon may 
have) appeared to be poorly described by the theory even when Sudakov resummations were carried out . This led 
£N) ' to some authors to claim that a non- vanishing gluon mass was necessary in order to describe the data [j| , even when 
relativistic corrections were taken into account Later on, with the advent of Non-Relativistic QCD (NRQCD) Q, 
these decays could be analyzed in a framework where short distance effects, at the scale of the heavy quark mass m 
or larger, could be separated in a systematic manner These short distance effects are calculated perturbatively 
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in a s (m) and encoded in matching coefficients whereas long distance effects are parameterized by matrix elements of 



local NRQCD operators. Even within this framework, a finite gluon mass seemed to be necessary to describe data 
l/"} , 0- However, about the same time it was pointed out that in the upper end-point region the NRQCD factorization 
approach breaks down and shape functions, namely matrix elements of non-local operators, rather than NRQCD 

^3 



matrix elements, must be introduced Q^l . Early attempts to modeling color octet shape functions produced results 
in complete disagreement with data [ljj, and hence later authors did not include them in their phenomenological 
I ; analysis [Hill. Notwithstanding this region has received considerable attention lately, as it was recognized that the 
so called Soft-Collinear Effective Theo ry (SC ET)L4, Il5 | may help in organizing the calculation and in performing 
resummations of large (Sudakov) logs I n fact, the early resummation of Sudakov logarithms Q has 



been recently correct ed I l7l within this framework, and statements about the absence of Sudakov suppression in the 
color singlet channel |18| have been clarified [13 ■ For the T(15) state, the bound state dynamics is amenable of a 
weak coupling analysis, at least as far as the soft scale (tuv. v ~ a s {mv) « 1, the typical velocity of the heavy quark 
in the quarkonium rest frame) is concerned pH I20H21II22I I23L l2^ | . These calculations can most conveniently be done 
in the framework of potential NRQCD (pNRQCD), a, further effective theory where the contributions due to the soft 
and ultrasoft (~ mv 2 ) scales are factorized |25ll2^.l27| (see [^3 for a review). Recently a calculation of the color octet 
shape functions, which combines SCET and pNRQCD, has become available [29) . 

It is the aim of this work to put together all known theoretical ingredients for these decays in order to see if a good 
description of data is achieved in the whole range of z, without the introduction of a finite gluon mass [30j |. The 
theoretical calculation of the so called direct contributions is under good parametric control in the central region and 
in most of the upper end-point region. Indeed, if one uses the original NRQCD velocity counting rules Q, namely 
a s (m) ~ v 2 and a s (mv) ~ v, together with existing calculations at weak coupling, a complete NLO expression can 
be put forward in the central region. For the upper end-point region a complete LO expression, which includes color 
octet contributions and takes into account both Sudakov and Coulomb resummations, is also available. The merging 
of the central and the upper end-point regions will be discussed in detail. The fragmentation contributions, i. e. those 
for which the photon originates from the decay products of the heavy quarks, are parametrically of the same order 
as the direct photon contributions in the central region [3lT ] and overweight the direct photon contributions in the 
lower end-point region {z — * 0). They play a minor, but non- negligible, role in our analysis and are subject to large 
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theoretical uncertainties. 

We distribute the paper as follows. In the next section we separate the contributions to the decay width into direct 
and fragmentation. Sections IIIII and llVI are devoted to either contributions respectively. In section we carry out 
the phenomenological analysis and section lv"fl is devoted to the conclusions. 



II. THE PHOTON SPECTRUM 

The contributions to the decay width can be split into direct ( dlr ) and fragmentation (f ra 9) 

dT _ dT dlr + dT /ra9 
dz dz dz 

We will call direct contributions to those in which the observed photon is emitted from the heavy quarks and fragmen- 
tation contributions to those in which it is emitted from the decay products (light quarks). This splitting is correct 
at the order we are working but should be refined at higher orders, z E [0, 1] is defined as z — 2E 1 /M (M is the 
mass of the heavy quarkonium state) , namely the fraction of the maximum energy the photon may have in the heavy 
quarkonium rest frame. 



III. DIRECT CONTRIBUTIONS 

The starting point is the QCD formula 

tjr dir M f 

^— = z— ImT(z) T{z) = -iJ foe-*** (V Q (nS)\T{J,(x)M0)}\V Q (nS)) rtf (2) 

where J^{x) is the electromagnetic current for heavy quarks in QCD and we have restricted ourselves to 3 S± states. 
q is the photon momentum, which in the rest frame of the heavy quarkonium is q = (q+ , q- , q± ) = [zM, 0,0). We 
have used light cone coordinates q± — q° ± q 3 . The approximations required to calculate J5J) are different in the lower 
end-point region {z — ► 0), in the central region (z ~ 0.5) and in the upper end-point region {z — > 1). We will denote 
Y dlr by r c and r e in the central and upper end-point regions respectively (the expressions for the lower end-point 
region will not be necessary). 



A. The central region 

For z away from the lower and upper end-points (0 and 1 respectively), no further scale is introduced beyond those 
inherent of the non-relativistic system. The integration of the scale m in the time ordered product of currents in 
leads to local NRQCD operators with matching coefficients which depend on m and z. At leading order one obtains 

1 dVfr, 2-z z(l-z) 1-z , s (1-z) 2 / 

rwf = — + (hw + 2 — ln(1 " z) " 2 W^f ln(1 " z) ' (3) 

where 

T = -aa s e Q , (4) 

and eQ is the charge of the heavy quark. The a s correction to this rate was calculated numerically in ref. |32l |. The 
expression corresponding to (0J in pNRQCD is obtained at lowest order in any of the possible regimes by just making 
the substitution 

(VQinS^OxCSi^VQinS)) = 27V c |^„ (O)| 2 , (5) 

where V'no(O) is the wave function at the origin. The final result coincides with the one of the early QCD calculations 
0,0. We will take the Coulomb form Vio(O) = 7 3 /tt for the LO analysis of T(15) (7 is defined in l|AT3)l ). 
The NLO contribution in the original NRQCD counting is v 2 suppressed with respect to (j2J- It reads 
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^ LO = (3 v (V Q (nS)\V 1 ( 3 S 1 )\V Q (nS)) 
dz l[ - 1} m 4 1 j 

In the original NRQCD counting or in the weak coupling regime of pNRQCD the new matrix element above can 
be written in terms of the original one |33| 1 

(V Q (nS)\V i eS 1 )\V Q (nS)) = / M - 2m \ (j^(ng)jOi(^)jj^(n^ ^ | Q , 2)) 
m 4 \ m J m 2 ^ \ JJ 

The matching coefficient can be extracted from an early calculation [(| (see also [35|L It reads 

C t ( 3 Si) = -^aa 2 s e 2 Q (f b (z) + l -F w (z)^ (8) 

where (£ = 1 — z) 

2 - 16^ + 10<f - 48£ 3 - 10£ 4 + 64£ 5 - 2£ 6 + (1 - 3£ + 14£ 2 - 106<f + 17£ 4 - 5Lf) ln£ 



Fb(z) = 



F w (z) 



2(l-0 3 (l + 4 

-26 + 14£ - 210C 2 + 134£ 3 + 274£ 4 - 150f 5 - 38£ 6 + 2£ 7 
= 3(l-C) 3 (l + 5 

(27 + 50£ + 257£ 2 - 292£ 3 + 205£ 4 - 78£ 5 - 41£ 6 ) ln£ 
3(l-0 3 (l + 5 



(9) 



The contributions of color octet operators start at order v 4 . Furthermore, away of the upper end-point region, 
the lowest order color octet contribution identically vanishes ||. Hence there is no l/a s enhancement in the central 
region and we can safely neglect these contributions here. 

If we use the counting a s (fih) ~ v 2 , a s (/i s ) ~ v (nh ~ m and fi s ~ mv are the hard and the soft scales respectively) 
for the T(15), the complete result up to NLO (including v 2 suppressed contributions) can be written as 

= dT|o + dT c NLO + dVlp^ 
dz dz dz dz 

The first term consist of the expression Q with the Coulomb wave function at the origin (JSJ including corrections 



up to O (a s (fis)) [36L l37j|. the second term is given in 10, and the third term consists of the radiative O (a s (nh)) 

corrections to (jSJ which have been calculated numerically in |32| . Let us mention at this point that the O (a s (/i s )) 2 

corrections to the wave function at the origin turn out to be as large as the leading order term. This will be important 
for our final results. Note that the standard NRQCD counting we use does not coincide with the usual counting of 
pNRQCD in weak coupling calculations, where a s (/x/ t ) ~ a s (/i s ) ~ a s (mv 2 ). The latter is necessary in order to get 
factorization scale independent results beyond NNLO for the spectrum and beyond NLO for creation and annihilation 
currents. However, for the T(1S) system (and the remaining heavy quarkonium states) the ultrasoft scale mv 2 
is rather low, which suggests that perturbation theory should better be avoided at this scale |2(j. This leads us 
to standard NRQCD counting. The factorization scale dep endences that this counting induces can in principle be 
avoided using rcnormalization group techniques 3.8, 39|. |40l l4ll [42j . In practice, however, only partial NNLL results 
exists for the creation and annihilation currents [Ij, [lij (see j45j for the complete NLL results), which would fix 
the scale dependence of the wave function at the origin at 0(a 2 (mv)). We will not use them and will just set the 
factorization scale to m. 



1 In the strong coupling regime of pNRQCD an additional contribution appears |34 
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B. The lower end-point region 



For z — > 0, the emitted low energy photon can only produce transitions within the non-relativistic bound state 
without destroying it. Hence the direct low energy photon emission takes place in two steps: (i) the photon is emitted 
(dominantly by dipole electric and magnetic transitions) and (ii) the remaining (off-shell) bound state is annihilated 
into light hadrons. It has a suppression ~ z 3 with respect to Tq (see f° r a recent analysis of this region 

in QED). Hence, at some point the direct photon emission is overtaken by the fragmentation contributions l3l|. 
In practice this happens about z ~ 0.4, namely much before than the z 3 behavior of the low energy direct photon 
emission can be observed, and hence we shall neglect the latter in the following. 



C. The upper end-point region 

In this region the standard NRQCD factorization is not applicable This is due to the fact that small scales 
induced by the kinematics enter the problem and have an interplay with the bound state dynamics. In order to study 
this region, one has to take into account collinear degrees of freedom in addition to those of NRQCD. This can be 
done using SCET as it has been described in [13l llq . In this region, the color octet contributions are only suppressed 
by v 2 or by 1 — z. Since their matching coefficients are enhanced by l/a s (^h), they become as important as the color 
singlet contributions if we count a s (nh) ~ v 2 ~ 1 — z. We will write 

dT e dT e cs dT e co 

= ^£ + _2£ (11) 

az az az 

where CS and CO stand for color singlet and color octet contributions respectively. 

1, Color singlet contributions 
For the color singlet contribution we shall use the expression with the Sudakov resummed coefficient in ref. 



Z n odd [ J5/2 



+ 



3/ 



(n) 



3/2 Igq 



(12) 



where the definitions for the different functions appearing in i|12|) are collected in the Appendix 1X1 



2. Color octet contributions 



For the color octet contributions we use 

M 

dT co , , , N /l6Ma\ . NJ 

dz = a s {Hu)a s (nh) I 81to4 J J C(x - z)b s +p{x)dx (13) 

/j, u is the ultrasoft (US) scale (see section^for the expression we use). C(x — z) contains the Sudakov resummations 
of ref. [3 2 , 

r( _ ] = _±L ( _ ) exp[^ 1 [a s ^/(4 7 r)]+g 2 [a s /3 £/(4 7 r)]] 1 

1 ' dz\ [ ' r[l-« ?1 [a s /3o^/(47r)]~a s/ 3o^/(47r) 5 ' 1 [a s /3n^/(4^)]]J 1 ' 



2 The matching coefficients provided in this reference become imaginary for extremely small values of z — 1, a region where our results do 
not hold anyway. We have just cut-off this region in the convolutions. 
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The (tree level) matching coefficients (up to a global factor) and the various shape functions are encoded in Ss+p(x), 

•<«> : = " {- {^WT (S) 2 P S(M(1 - >» - {^wY 1 <3S " (J,(1 - + - •)») < 15 > 



>S+P{ 



The definitions of all the functions appearing in l|14fl and (|15ll are collected in the Appendix [S] The shape functions 
Ss, Spi and S P2 may become Sg IS , Sp{ s and Sp 2 s or Sg ub , Spf and Spf depending on the subtraction scheme 
employed. The procedure used to renormalize the shape functions is explained in the Appendix^] In figure^we plot 
the end-point contribution with the shape functions renormalized in an MS scheme (blue dashed line) and in the 
sub scheme, which makes additional subtractions (red solid line), together with the experimental data |48| (we have 
convoluted the theoretical curves with the experimental efficiency, the overall normalization of each curve is taken as 
a free parameter. For the details of the scale setting see section 01. We see that both schemes are equally good for 
the description of the shape of the experimental data in the end-point region. 

Note (from Appendix A) that we use the octet shape functions calculated in ref. [23, which are crucial in order 
to have a good description of data in the upper end-point region. We would like to comment on the validity of those 
formulas. This is limited by the perturbative treatment of the US gluons. The typical momentum of these gluons in 
light cone coordinates turns out to be: 



(k+,k u fc_) = I M(l - z), J2M(1 - z) [ M{ \ Z) ~EA, M(l - z) - 2Bh I ( U> ) 



Note that the typical k± is not fixed by the bound state dynamics only but by a combination of the latter and the 
end-point kinematics. Hence, the calculation is reliable provided that k± > lGeV., which means z < 0.92. Note also 
that the typical three momentum of the heavy quarks in the shape function is given by 



P~^m(j(l-z)-Eh) (17) 

This means that the multipole expansion always holds in the end-point re gion , regardless that M(l — z) is bigger or 
smaller than mal . This important point was not sufficiently emphasized in [29| and it is the ultimate reason why such 
a good description of data was obtained there for z € [0.7, 1]. Recall that at z < 0.7 or so we are at the border of the 
end-point region and the standard NRQCD factorization formulas should hold. This means that the contributions of 
the octet shape functions should merge suitable contributions of the NRQCD factorization formulas, as we discuss 
below. 

D. Merging the central and upper end-point regions 

As we have seen, different approximations are necessary in the central and upper end-point regions. It is then not 
obvious how the results for the central and for the upper end-point regions must be combined in order to get a reliable 
description of the whole spectrum. When the results of the central region are used in the upper end-point region, 
one misses certain Sudakov and Coulomb resummations which are necessary because the softer scales My/1 — z and 
M(l — z) become relevant. Conversely, when results for the end-point region are used in the central region, one misses 
non-trivial functions of z, which are approximated by their end-point [z ~ 1) behavior. 

One way to proceed is the following. If we assume that the expressions for the end-point contain the ones of the 
central region up to a certain order in (1 — z), we could just subtract from the expressions in the central region the 
behavior when z — > 1 at the desired order and add the expressions in the end-point region. Indeed, when z -t 1 
this procedure would improve on the central region expressions up to a given order in (1 — z), and when z belongs 
to the central re gion , they would reduce to the central region expressions up to higher orders in a s . This method 
was used in ref. |l3j ] and in ref. j4^|. In ref. only color singlet contributions were considered and the end-point 
expressions trivially contained the central region expressions in the limit z — » 1. In ref. j4fj color octet contributions 
were included, which contain terms proportional to (1 — z). Hence, the following formula was used 

1 dT dlr 1 dT L jn / 1 dT e nq \ / 1 dT e rn \ 

^ ->- 1 ^ - z ) + [ -fQ- - z (4 + 2 log (1 - z)) (1 - z)\ (18) 



ro dz To dz \To dz I \ To dz 



Even though a remarkable description of data was achieved with this formula (upon using a suitable subtraction 
scheme described below), this method suffers from the following shortcoming. The hypothesis that the expressions 
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for the end-point contain the ones for the central region up to a given order in (1 — z) is in general not fulfilled. As 
we will see below, typically, they only contain part of the expressions for the central region. This is due to the fact 
that some a s {pkh) in the central region may soften as a s (M(l — z)), others as a s (My/l — z) and others may stay at 
Qtsifi'h) when approaching the end-point region. In a LO approximation at the end-point region, only the terms with 
the a s at low scales would be kept and the rest neglected, producing the above mentioned mismatch. We shall not 
pursue this procedure any further. 

Let us look for an alternative. Recall first that the expressions we have obtained for the upper end-point region are 
non-trivial functions of M(l — z), My/1 — z, ma s (mv) and ma 2 (mv), which involve a s at all these scales. They take 
into account both Sudakov and Coulomb rcsummations. When z approaches the central region, we can expand them 
in a s (My/l — z), a s (M(l — z)) and the ratio ma s (mv)/My/l — z. They should reduce to the form of the expressions 
for the central region, since we are just undoing the Sudakov and (part of) the Coulomb resummations. Indeed, we 
obtain 



dT e C s 




dz 


dz 


dTho 


, aL CO 


dz 


dz 



T z (l + j± (CU (2vr 2 - 17) + 2n f ) log(l - z) + 0(a 2 s ) 
16Ma\ 



81m 4 

+M(l-z) [-1 + log 



J 2 |-0io (0)| ( ma s VT^A + 



i4 



M 2 (l-zf 



-M- 



2tt 
23 



1 



-2C A [ ^(l-z)log 2 (l-z) 



|c A -^j((l-z)log(l-z) 



log 



log 



M 2 (l-z)- 

2 



t'c 



-^K l0S (w^F 



M 2 (l-zf 
7 2 7 4 

I) + [ ma sl a s — , -^r 
m m° 



(19) 



dx — — ^-f{x 1 z) 



1 



dx f(x, z) 

x — z 



(20) 



where 



f(x,z) = (1 - x)log 



M 2 (l ~xf 



- (1 - z) log 



M 2 (l - z)< 



(21) 
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A = —N c — 136C/(2 — A)/9 (in an MS scheme; it becomes A = — 64C/(2 — A)/9 in the sub scheme described in 
Appendix B). The details of this derivation are given in the Appendix [CJ The color singlet contribution reproduces 
the full LO expression for the central region in the limit z — > 1. The color octet shape functions Spi and Sp2 give 
contributions to the relativistic corrections JHJ), and Sp2 to terms proportional to (1 — z) in the limit z — > 1 of @ as 
well. We have checked that, in the z — * 1 limit, both the (1 — z) ln(l — z) of © and the ln(l — z) of the relativistic 
correction © are correctly reproduced if [i c ~ My/1 — z, as it should. All the color octet shape functions contribute 
to the 0(a s (nh)) correction in the first line of (|2(J|) . There are additional 0(a s (/j,h)) contributions coming from the 
expansion of the (Sudakov) resummed matching coefficients of the color singlet contribution and of the Sp2 color 
octet shape function. The a s log(l — z) in (|T9"|) reproduces the logarithm in dT c Lo,a B /dz. 
We propose the following formula 



1 dT di > 
Tq dz 



i dr c 

Tq dz 



1 dT 



cs 



1 dT 



cs 



Tq dz 



dz 



1 dT 



CO 



1 dT 



CO 



Tq dz Tq dz 



(22) 



This formula reduces to the NRQCD expression in the central region. When we approach the upper end-point region 
the second terms in each of the parentheses are expected to cancel corresponding terms in the z — » 1 limit of the 
expression for the central region up to higher order terms (in the end-point region counting). Thus, we are left with 
the resummed expressions for the end-point (up to higher order terms). 

There are of course other possibilities for the merging. For instance, one may choose a z\ below which one trusts 
the calculation for the central region and a z-i above which one trusts the end-point region calculation, and use some 
sort of interpolation between z\ and Z2 (see for instance [50]). This would have the advantage of keeping the right 
approximation below z\ and beyond z 2 unpolluted, at the expense of introducing further theoretical ambiguities due 
to the choice of z\ and Z2, and, more important, due to the choice of the interpolation between z\ and zi- We believe 
that our formula ll'L'l) is superior because it does not introduce the above mentioned theoretical ambiguities. The price 
to be paid is that the expressions from the central region have an influence in the end-point region and vice- versa. This 
influence can always be chosen to be parametrically sublcading but large numerical factors may make it noticeable in 
some cases, as we shall see below. 



1. Merging at LO 

If we wish to use only the LO expressions for the central region, we should take l(20jl at LO, namely 

I dT 



1 dT* cs 



Tq dz 



CO 



Tq dz 



and substitute them in l|22|) . Unexpectedly, the results obtained with this formula in the central region deviate 
considerably from those obtained with formula © (see Fig. [2J|. This can be traced back to the fact that the a s y/l — z 
corrections in l|2l)|l are enhanced by large numerical factors, which indicates that the merging should better be done 
including a s (fih) corrections in the central region, as we discuss in the next section. Alternatively, we may change our 
subtraction scheme in order to (partially) get rid of these contributions, as discussed in the Appendix El With the 
new subtraction scheme (sub) the situation improves, although it does not become fully satisfactory (see Fig. |2J). This 
is due to the fact that some a s y/l — z terms remain, which do not seem to be associated to the freedom of choosing a 
particular subtraction scheme. In spite of this the description of data turns out to be extremely good. In figure 01 we 
plot, using the sub scheme, the merging at LO (solid red line) and also, for comparison, equation l|18|l (blue dashed 
line). We have convoluted the theoretical curves with the experimental efficiency and the overall normalization is 
taken as a free parameter. 



2. Merging at NLO 

If we wish to use the NLO expressions for the central region (|I0|l , we should take all the terms displayed in l|19fl- 
(120(1 and substitute them in l(22|l . Unlike in the LO case, for values of z in the central region the curve obtained from 
(|22|l now approaches smoothly the expressions for the central region I|1U|) as it should. This is so no matter if we 
include the a^(fj. s ) corrections to the wave function at the origin in dT c LO /dz, as we in principle should, or not (see 
Figs. [Hand|SJ. However, since the above corrections are very large, the behavior of the curve for z — ► 1, strongly 
depends on whether we include them or not (see again Figs. 0] and [5J . We believe that the two possibilities are 
legitimate. If one interpretes the large a^(u s ) corrections as a sign that the asymptotic series starts exploding, one 
should better stay at LO (or including a B (/x s ) corrections). However, if one believes that the large a^(fj, s ) corrections 
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FIG. 2: Merging at LO. The solid red line is the NRQCD expression The dot-dashed curves are obtained using an MS 
scheme: the pink (light) curve is the end-point contribution (li lt and the black (dark) curve is the LO merging of section llll D II 
The dashed curves are obtained using the sub scheme (explained in the Appendix [EJ: the green (light) curve is the end-point 
contribution 1111 and the blue (dark) curve is the LO merging of section fill D II 




FIG. 3: Direct contribution to the spectrum. The solid red line corresponds to the LO merging of section fill D II and the blue 
dashed line corresponds to equation 1181 . The points are the CLEO data |48|. 
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FIG. 4: Merging at NLO (using an MS scheme and the wave function at the origin at LO). The solid red line is the NRQCD 
result (I1UL the blue (light) dashed curve is the end-point contribution 11 U and the black (dark) dashed curve is the NLO 
merging of section lill D 21 



are an accident and that the a^(/j, s ) ones (see |5ll 15^ for partial results) will again be small, one should use these 
ckg(/i s ) corrections. We consider below the two cases. 

If we stay at LO (or including a s (p s ) corrections) for the wave function at the origin, the curve we obtain for 
z — *■ 1 differs considerably from the expressions for the end-point region (see Fig. 0J. This can be traced back 
to the a s \/\ — z term in (|20J) again. This term is parametrically suppressed in the end-point region, but, since it 
is multiplied by a large numerical factor, its contribution turns out to be overwhelming. This term might (largely) 
cancel out against higher order contributions in the end-point region, in particular against certain parts of the NLO 
expressions for the color singlet contributions, which are unknown at the moment. 

If we use the wave function at the origin with the a!?(/i s ) corrections included, the curves we obtain for z — > 1 
become much closer to the expressions for the end-point region Ijlll) (see Fig[SJ). Hence, a good description of data 
is obtained with no need of additional subtractions, as shown in figure |H| (as usual experimental efficiency has been 
taken into account and the overall normalization is a free parameter). 



IV. FRAGMENTATION CONTRIBUTIONS 



The fragmentation contributions can be written as 

dT frag nl fa 



dz 



£ / -C a (x)D a J*,M), (24) 



where C a represents the partonic kernels and D ai represents the fragmentation functions. The partonic kernels can 
again be expanded in powers of v [j| 

C a =J2 C aiQ] (25) 
Q 

The leading order term in v is the color singlet rate to produce three gluons 

c, [0l e a)] - ^ + - .) - J£*hi - .)) (26) 
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FIG. 5: Merging at NLO (using an MS scheme and the wave function at the origin with the ag(p s ) corrections included). The 
solid red line is the NRQCD result dlUt . the blue (light) dashed curve is the end-point contribution (111! and the black (dark) 
dashed curve is the NLO merging of section lrTTTT2l 




FIG. 6: Direct contribution to the spectrum using the NLO merging of section lTlID2I C in an MS scheme and the wave function 
at the origin with the ag(/x s ) corrections included). The points are the CLEO data 48|. 
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The color octet contributions start at order u 4 but have a — enhancement with respect to (I26|) 

rlnno.i ^a 2 ^ (y Q (n5)|0 8 (^ o )|y Q (n5)) 
C g [U s [ S )\ - — — d(l - z) —2 

3 m 4 

C q [O s (* Sl) ] = 5^5(1 - Z) <M«M (27) 

3 

Then the color singlet fragmentation contribution is of order a 3 D ff _> 7 and the color octet fragmentation are of 
order v A a 2 Dg^ ( 1 S , and 3 Pj contributions) or v A a 2 s D q _~ f ( 3 S\ contribution). We can use, as before, the counting 
v 2 ~ a s to compare the relative importance of the different contributions together with the existing models for the 
fragmentation functions |53j | . The latter tell us that Z? q ^ 7 is much larger than Z? g ^ 7 . This causes the O^a^Dq^) 
3 Si octet contribution to dominate in front of the singlet 0(a 3 D g ^ 1 ) and the octet O^a^Dg^y) contributions. 
In fact, a s D q ^^ is still larger than D 9 _> 7 , so we will include in our plots the a s corrections to the color octet 
contributions (|27[1 proportional to D q ^, which have been calculated in 8]. In addition, the coefficients for the 
octet 3 Pj contributions have large numerical factors, causing these terms to be more important than the color singlet 
contributions. Let us finally notice that the a s corrections to the singlet rate will produce terms of C(a 4 Z? 9 _, 7 ), which 
from the considerations above are expected to be as important as the octet 3 Si contribution. These a s corrections to 
the singlet rate are unknown, which results in a large theoretical uncertainty in the fragmentation contributions 

For the quark fragmentation function we will use the LEP measurement 



D ql (z,fi) 



e|o(M) 
2tt 



(28) 



where 



M 2 1 4- fl — 7^ 
C = -1-1h(tt) ; *Vr(*) = — L ; Mo = 0.14t°;g GeV (29) 



and for the gluon fragmentation function the model [F^. These are the same choices as in JT|j. However, for the 
^8( 1 5'o) an d 0$( 3 Pq) matrix elements we will use our estimates in [2^ 

<T(15)|O 8 ( 1 5 )|T(15))| M=M ~ 0.004Gey 3 (30) 
(r(lS)|0 8 ( 3 Po)|T(lS))| M=M ~ 0.08 GeV 5 (31) 

The above numbers are obtained in an MS scheme from dimensionally regularized US loops only. The value 
we assign to the S'-wave matrix element is compatible with the recent (quenched) lattice determination (hybrid 
algorith m') [ 56l . Notice that we do not assume that a suitable combination of these matrix elements is small, as it was 
done in |l3j|. The Og,( 3 Si) matrix element can be extracted from a lattice determination of the reference above |56| . 
Using the wave function at the origin with the ag(/j, s ) corrections included, we obtain, 

(T(1S*)|C 8 ( 3 5i)|T(1S'))| ai=a/ - 0.00026GeU 3 (32) 

which differs from the estimate using NRQCD v scaling by more than two orders of magnitude: 

(T(15)|0 8 ( 3 5 1 )|T(15))|^ M ~ v 4 <T(15)|0 1 ( 3 5 1 )|T(15))| ai=m ~ 0.02 GeV 3 (33) 

(we have taken v 2 ~ 0.08), which was used in ref. 0. The description of data turns out to be better with the 
estimate l|33ll . However, this is not very significant, since, as mentioned before, unknown NLO contributions are 
expected to be sizable. 

In the z — > region soft radiation becomes dominant and the fragmentation contributions completely dominate 
the spectrum in contrast with the direct contributions |3l| . Note that, since the fragmentation contributions have an 
associated bremsstrahlung spectrum, they can not be safely integrated down to z = 0; that is Jq dz dr ^ z 9 is not an 
infrared safe observable. In any case we are not interested in regularizing such divergence because the resolution of 
the detector works as a physical cut-off. 
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V. SCALE SETTING AND ERROR ANALYSIS 



Formula 1)221) requires dT e /dz for all values of z. The color octet shape functions, however, were calculated in the 
end-point region under the assumption that My/1 — z ~ 7, and the scale of the a s was set accordingly. When z 
approaches the central region My/1 — z 7, and hence some a s will depend on the scale My/1 — z and others on 
7 (we leave aside the global a s (/•*«)) which will be discussed below). In order to decide the scale we set for each a s 
let us have a closer look at the formula (|20[1 . We see that all terms have a common factor 7 3 . This indicates that 
one should extract 7 3 factors in the shape functions, the a s of which should stay at the scale /i s . T his is achieved by 
extracting 7 s / 2 in 1$ and Ip. If we set the remaining a s to the scale \i v = y/m(M{l — z)/2 — E{), we will reproduce 
(H3> when approaching to the central region, except for the relativistic correction, the a s of which will be at the scale 
Hp instead of at the right scale fi s . We correct for this by making the following substitution 

472 _ mCja 2 M \ 
3m 3 J 

Notice that the replacements above are irrelevant as far as the end-point region is concerned, but important for the 
shape functions to actually (numerically) approach the expressions (|20|) in the central region, as they should. 

We use the following values of the masses for the plots: m — 4.81 GeV and M = 9.46 GeV. The hard scale 
/i/j is set to fih = M. The soft scale fi s = mCfa s is to be used for the a s participating in the bound state 
dynamics, we have a s (fj, s ) = 0.28. The US scale fi u , arising from the couplings of the US gluons, is set to (i u — 

yj2M(l — z) (4f(l — z) — Ei) (as discussed in section ITTTf) . We have used the Mathematica package RunDec [57| to 

obtain the (one loop) values of a s at the different scales. 

Our final plot in Fig. 0is obtained by using the merging formula (|22|) at NLO with the a 2 (/Lt s ) corrections to the 
wave function at the origin included for the direct contributions plus the fragmentation contributions in section llVl 
including the first a s corrections in C q and using the estimate 1)33)1 for the (T(lS')|08( 3 5'i)|T(lS')) matrix element. 

The error band is obtained by replacing fi c by V2 ±:L /i c . Errors associated to the large a 2 (/i s ) corrections to the 
wave function at the origin, to possible large NLO color singlet contributions in the end-point region and to the 
fragmentation contributions are difficult to estimate and not displayed (see the corresponding sections in the text for 
discussions). The remaining error sources are negligible. As usual experimental efficiency has been taken into account 
and the overall normalization is a free parameter. 



VI. CONCLUSIONS 



We have analyzed the photon spectrum in radiative T(IS') decays within an Effective Field Theory framework. For 
the direct contributions, the merging of the results for the central and upper end-point regions has been discussed 
in detail. We have shown how to consistently combine the complete LO results for the upper end-point region with 
the complete NLO ones for the central region. We have seen that the large a 2 (/i s ) corrections to the wave function 
at the origin are important in order to get a good description of data. Otherwise, parametrically subleading large 
contributions in the end-point region would be necessary. We would like to emphasize that our final results for 
the direct contributions are essentially parameter free: only the mass of the bottom quark to, the strong coupling 
constant a s , and the proper choice of subtraction scales (which appear in logarithms) are used as an input. For the 
fragmentation contributions, we have pointed out that if the commonly used model for the gluon fragmentation into 
a photon is appropriated, a s corrections to the LO color singlet matrix element giving rise to a light quark which 
fragments into a photon may be as important as the LO results. Hence, fragmentation contributions suffer from large 
theoretical uncertainties. Nevertheless, if we put together the available theoretical results for these contributions with 
the ones for the direct contributions, an excellent description of data is achieved for the whole part of the spectrum 
where experimental errors are reasonable small. Clearly, our results indicate that the introduction of a finite gluon 
mass | 3(| is unnecessary. One should keep in mind, however, that in order to have the theoretical errors under control 
higher order calculations are necessary both in the direct (end-point) and fragmentation contributions. 

Before closing, let us mention that the inclusion of color octet contributions in the end-point region together with 
the merging with the central region expression described in this work may be useful for production processes like 
inclusive J ftp production in e + e~ machines |50l l58l 159^ . 
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FIG. 7: Photon spectrum. The points are the CLEO data |4ftj |. The solid lines are the NLO merging in section HHP 21 
plus the fragmentation contributions: the red (light) line and the blue (dark) line are obtained by using 1331 and 1321 for 
(T(1S')|08( 3 S'i)|T(1S')) respectively. The grey shaded region is obtained by varying fi c by y/2 ±1 fi c - The green shaded region 
on the right shows the zone where the calculation of the shape functions is not reliable (see section llll C 21 . The pink dashed 
line is the result in 0]> where only color singlet contributions were included in the direct contributions. 
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APPENDIX A: DEFINITIONS 

In this appendix we collect the definitions for the formulas that appear in the paper. 

C N C N *~ l R (ng.4 - 2n f ) 34Cj - 10C A n f - 6C f n f 
Ca ~ Nc C/ = ^VT 00 = 3 A = 3 

N c = 3 is the number of colors, rif = 4 is the number of light flavors. 



1. Definitions for formula (1121) 



(„) _ n(n + l)(n + 2)(n + 3) 



5/2 



9(n + 3/2) 



(n) _ (n + l)(n + 2) 



3/2 



■3/2 



(Al) 



"(/*) = 



a s (fJ-) 
a s (2m) 



(n) = 7gg ~ A T 
r± A 



(n) 
Iqq 



'99 



'99 



~(») 
'99 



= (7/ 



1 



i 11+1 1 
-i-2Vi 



(n + l)(n + 2) 2 
n 2 + 3n + 4 



3 C/ (r/ - I'll// + 2) 



3n/ 

cu 



3n + 4 



n(n + l)(ra + 2)(ra + 3) 
2 2 



1 



n(n+l) (n + 2)(n + 3) 6 



n+l 



9i (x) 
52 (x) 



2. Definitions for formula (11 111 



£ Ri — log(x — 2) 



2r 



adj 



^ [(1 - 2 X ) log(l - 2 X ) - 2(1 - x) log(l - x)] 



/3oX 



8r 2 



adj 



[-log(l-2x)+21og(l- X )] 



2Tf% 

PI 



log(l - 2 X ) - 21og(l - X) + g iog^ 1 - - W(l - X) 



+ log(l - X) + -t 1 log(l - 2x) - —j- log no [log(l - 2 X ) - log(l - X )] 

Po Po Po 



c A ; r 



CU 



67 7T 2 

36 ~ 12 



5n/ 
18 



B X = -C A ; 7l 



^0 
' 4 



n = e 



"7E 



3. Definitions for formula (115 [I 



,fc+ /t a s N c 1 2z' „ / A _ A 1 - z' 



7s(— +x) = m 



7T 2 1 - z' 



1 - 



M^T + *) = V V 3 (2 " A) 4(l + z0 3 i 2(1 + Z ' K2 + ^ + (5 + 3Z ' K_1 + A) + 2M + A 
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S ^ = W (B 2 f * ( 2 ^(0)/,(^ + x) - lU k f + x)) (MO) 

= ^ jf* 5 (Vo(O)M^ + x) - #(*± + *)) (All) 

S P2 (k + ) = f dx {k 8k +2 x f (^» " 2 ^o(0)M^ + x) + 4(^ + a:)) (A12) 

« .^^.^ A = --i— E 1 = -^ (A13) 
'2 7 m 2 27V C C/ m y ' 

cf is the hard matching coefficient of the chromomagnetic interaction in NRQCD, it will be taken to 1. The renor- 
malized expressions (see Appendix iBll in an MS scheme read 

- W<£)'{"-»(-/P?) (r( r * ( T-)4-<-»^ h ^)*- 



c 2 F a s (n u )"f 3 C 2 f a 2 s (n P ) , -nW , ( ^ 



3Tr 2 N c m 



(l-A)(-2 + A(21n2 + l)) In (A14) 



m ^- + ->- 

1 2 1 m 



6ttJV c I 7 w \ V 7T 3 V / \ -A) V 2 2z ' V 4 4 / z' : 

2 

7 



'<=+ 


f 7> _ 











(2 -A) 



a s (^)7«K(M ( 31 , A(41n2 + 19)_ A2(21n2+ M J /' 



9ir 2 N r V 6 v 6 ; v 6 



2a s (^«)7 5 ( ,( M 



9w 2 N c m V \k 



\ 2 ~ TO 

1" ( 7T I 1 (AW) 



S?/( , +) = ^)| A(0)fe+ (_ 2 + 21n (| 



^ (fc+ + 2a;) 2 ^-2^o(0)/p(^+*)+JjM^- + , i j } !AK») 
and in the scheme where additional subtractions are carried out 



snt+) . + ^ ( m 7i^) (A17) 



S{.?(* + ) = S«/(t + ) + ^2 ! M0) l/2if(2-A) -i t /t 

07TiV c \ V 7T 3 / \/TCh 



where 



S s P l b (k + ) = S A P ^(k + ) (A19) 

7 f k+ > / Ra s N c \ 1 fc+ 

/s( T +.x) := [rrhj-—) fc( T +*) 

7 P (^+x) := f^/?l (2_A) ) (A20) 
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APPENDIX B: REGULARIZATION AND RENORMALIZATION 

The shape functions (|A10(1 - I|A12() are ultraviolet (UV) divergent and require regularization and renormalization. In 
ref. |2^| it was pointed out that using dimensional regularization (DR) for the US loop only was enough to regulate 
them. In fact, the expressions HA10|) - (|A11|I implicitly assume that DR is used, otherwise linearly divergent terms 
proportional to V'io(O) would appear (which make (|A10|I - I|A11|I formally positive definite quantities). In addition an 
MS scheme was used to subtract the poles. Since it turns out that the final outcome strongly depends on the details 
of this subtraction, let us spell out the procedure carried out in ref. [2j|. In order to isolate the 1/e poles, Is and Ip 
were expanded up to 0(l/z' ) . The result was subtracted and added to the integrand of (|A10|) - (IA11() (for l|A12[l this 
is not necessary since the only divergent piece is independent of Ip). The subtracted part makes the shape functions 
finite. The added part contains linear and logarithmic UV divergencies. The 1/e (D = 4 — 2e) poles displayed in 
formulas (16) of ref. |2^, and eventually subtracted, were obtained by making dx — > dx(x//i)~ E in (|A10(l - i|A12(l . This 
was motivated by the fact that x ~ (kj_ being the transverse momentum of the US gluon) but differs from a 
standard MS scheme. Linear divergences are set to zero as usual in DR. 

We have used here a regularization and renormalization scheme which is closer to the standard one in pNRQCD 
calculations. We have regulated both the US loop and the potential loops (entering in the bound state dynamics) 
in DR. We have identified US divergencies by taking the limit D — > 4 in the US loops while leaving the potential 
loops in D dimensions [6(|- Potential divergencies are identified by taken D — > 4 in the potential loops once the US 
divergencies have been subtracted. It turns out that all divergencies in Sp2 are US and all divergencies in Ss are 
potential. Spi contains both US and potential divergencies. The potential divergences related with the bound state 
dynamics can be isolated using the methods of ref. |6lj |. The formulas corresponding to (16) of ref. |29l ] in this scheme 
read 

Cpa s (u u )7 3 C?a^(u„) 




2a ^"( 1 In )+...) (Bl) 



9ir 2 N c m \ e 



^(M,^?4^(± + 1 »(f )+■■•) (B2, 

For simplicity, we have set D = 4 everywhere except in the momentum integrals. fi p is defined in section /i c and 
fip C are the subtraction points of the US and potential divergencies respectively. If we subtract the 1/e poles and set 
/i c = M \/l — z and fi pc = y/m/i c we obtain exactly the same result as in ref. |29| for what the potential divergences 

is concerned 3 . For the US divergences there is a factor ln (^q-) of difference with respect to the previous scheme. 

In ref. an additional subtraction related to linear divergencies was made. This subtraction was necessary in 
order to merge smoothly with the results in the central region. We will also need this subtraction here when merging 
at LO, as discussed in subsection IIII Dl We use 




, which differ from the MS scheme by the subtraction of the second term in the square brackets. 



3 We assume that the correlation of scales advocated in l3Sl (see |42| for the implementation in our framework) must also be taken into 
account here. 
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APPENDIX C: THE SHAPE FUNCTIONS IN THE CENTRAL REGION 



When z approaches the central region from the upper end-point, the shape functions should reduce to matrix 
elements of NRQCD operators multiplied by the corresponding matching coefficients. We will see here that this is 
indeed the case. 

Let us first consider the S-wave octet shape function as defined in [2{| 



+*):=/ ^ 10 (x) 1 - J! fc + (CI) 

h a = p 2 /m + V , V a = a 8 /(2iV c |r|). When z approaches the central region, k + ~ M(l — z) ^> — E\ and the larger 
three momentum scale is Myl — z S> 7, the typical three momentum in the bound state. Therefore we can treat the 
Coulomb potential in i|Cl|) as a perturbation when it is dominated by this scale. It is convenient to proceed in two 
steps. First we write h = h s + (V — V s ), where h s — p 2 /m + V s , V s — — a s C//|r|, and expand V Q — V s . This allows 
to set h s — Ei to zero in the left-most propagator and makes explicit the cancellation between the first term in the 
series and the first term in (|C1|) . It also makes explicit that the leading term will be proportional to a s (M^/l — z). 
Second, we expand V s in h s = p 2 /m + V S . In addition, since M^/\ — z ^> 7, the wave function can be expanded 
about the origin. Only the first term in both expansion is relevant in order to get (|20|l . 
Consider next the P-wave shape functions as defined in j2£* 



,k+ , 1 



Mf + „ „ __ j mm ^ _ K _^ + x ) V) <C2) 



In order to proceed analogously to the S'-wave case, we have first to move the x 1 away from the wave function 
/p(^+aO=^io(0) + ^^ |d 3 x^ 10 (x)|^— — 



fe_L 1 

-T+ x 



2V l 



Ei + ^+x\ m J h D -Ei + ^- 



(C3) 



For the left-most propagators we can now proceed as before, namely expanding V — V s . Note that the leading 
contribution in this expansion of the second term above exactly cancels against the first term. Of the remaining 
contributions of the second term only the next-to- leading one (0(a s )) is relevant to obtain H2U|) . Consider next the 
leading order contribution in this expansion of the last term. It reads 




— k — viy o — — k — W4 ( c4 ) 

h -Ei + ^-+x h -Ei + ^+x) J 

Now we proceed as before with the left- most propagators, namely expanding V a — V s . The leading order contribution 
of the first term above produces the relativistic correction 0(v 2 ) of (|20|l . The next-to-leading contribution of this 
term and the leading order one of the second term are 0(a s ) and also relevant to 120fl . The next-to-leading order 
contribution of the last term in l|C3fl in the V — V s expansion of the left-most propagator is also 0(a s ) and relevant 
to flUJ. 
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